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Abstract 

Transport of the Brownian particles driven by Levy flights coexisting with subdiffusion in asym- 
metric periodic potentials is investigated in the absence of any external driving forces. Using the 
Langevin-type dynamics with subordination techniques, we obtain the group velocity which can 
measure the transport. It is found that the group velocity increases monotonically with the subd- 
iffusive index and there exists an optimal value of the Levy index at which the group velocity takes 
its maximal value. There is a threshold value of the subdiffusive index below which the ratchet 
effects will disappear. The nonthermal character of the Levy flights and the asymmetry of the 
potential are necessary to obtain the directed transport. Some peculiar phenomena induced by the 
competition between Levy flights and subdiffusion are also observed. The pseudo-normal diffusion 
will appear on the level of the median. 
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I. INTRODUCTION 



Directed Brownian motion induced by zero-mean non-equilibrium fluctuations in the ab- 
sence of macroscopic forces and potential gradients is presently under intense investigation^. 
This comes from the desire to understand unidirectional transport in biological systems^, 
as well as their potential technological applications ranging from classical non-equilibrium 
models 3 to quantum systems 4 -. A ratchet system is generally defined as a system that is 
able to transport particles in a periodic structure with nonzero macroscopic velocity in the 
absence of macroscopic force on average. Broadly speaking, ratchet devices fall into three 
categories depending on how the applied perturbation couples to the substrate asymmetry: 
rocking ratchets^, flashing ratchets 6 -, and correlation ratchets^. Additionally, entropic ratch- 
ets, in which Brownian particles move in a confined structure, instead of a potential, were 
also extensively investigated 8 . These studies on directed transport of the Brownian particles 
focused on the normal diffusion. 

However, anomalous diffusion has attracted growing attention, being observed in vari- 
ous fields of physics and related sciences^. Description of physical models in terms of Levy 
flights and subdiffusion becomes more and more popular— ~—. In the complex systems the 
distinct class of subdiffusion processes was reported in condensed phases*^, ecology^, and 
biology^. Superdiffusion driven by Levy flights is actually observed in various real systems 
and is used to model a variety of processes such as bulk mediated surface diffusion 13 , exciton 
and charge transport in polymers under conformational motion 14 , transport in micelle sys- 
tems or heterogeneous rocks^two-dimensional rotating flow 16 , and many others^. Goychuk 
and coworkers^ studied the subdiffusive transport in tilted periodic potentials and estab- 
lished a universal scaling relation for diffusive transport. Dybiec and coworkers^ studied 
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the minimal setup for a Levy ratchet and found that due to the nonthermal character of 
the Levy noise, the net current can be obtained even in the absence of whatever additional 
time-dependent forces. Del-Castillo-Negrete and coworkers^ also found the similar results 
in constant force-driven Levy ratchet. Rosa and Beims^ studied the optimal transport and 
its relation to superdiffusive transport and Levy walks for Brownian particles in ratchet 
potential in the presence of modulated environment and external oscillating forces. We also 
studied the transport of Brownian particles in the presence of ac-driving forces and Levy 
flights and multiple current reversals were observed 2 ^- 

Recently, much attention has been devoted to the competition between subdiffusion and 
Levy flights. The competition is conveniently described by the fractional Fokker-Planck 
equation with temporal and spatial fractional derivatives 9 . It is very difficult to see this 
competition in the framework of the fractional Fokker-Planck dynamics. Magdziarz and 
coworkers 2 ^ proposed a equivalent approach based on the subordinated Langevin method 
to visualize the competition on the level of sample paths as well as on the level of prob- 
ability density functions. Based on this approach Dybiec and coworkers 23 found that due 
to the competition between Levy flights and subdiffusion, the standard measure used to 
discriminate between anomalous and normal behavior cannot be applied straightforwardly. 
Koren and coworkers 2 ^ have investigated the first passage times in one-dimensional system 
displaying a competition between subdiffusion and Levy flights and found some peculiar 
phenomena. 

What will happen when the particles move in a ratchet potential subjected to subdiffusion 
and Levy flights? In order to answer this question we use the subordinated Langevin method 
proposed by the Magdziarz and coworkers 22 to investigate this competition in a minimal Levy 
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ratchet without any external driving forces. We emphasize on visualizing the competition 
on the level of the group velocity and diffusion and finding how this competition affects the 
directed transport. 



II. MODEL AND METHODS 

We consider the transport of the Brownian particles driven by Levy flights and subdiffu- 
sion in the absence of whatever additional time-dependent forces. The competition between 
Levy flights and subdiffusion in a ratchet potential V(x) can be described by the fractional 
Fokker-Planck equation with temporal and spatial fractional derivatives^ 1 ^ 



dp(x, t) 



D}~ a 



d V (x) D 



p(x,t), (1) 



dx 1] d\x\ 

where p(x, t) is the probability density for particles at position x and time t. The prime 
stands for the derivative with respect to the space variable x. D is the anomalous diffusion 
coefficient which describes the noise intensity in the subordinated process, rj denotes the 
generalized friction constant. Here oD t ~ a is the fractional of the Riemann-Liouville operator 
(0 < a < 1) defined through^ 



oD\- a g(t) = —- / {t - sr- 1 g(s)ds, (2) 



1 d /■*. 
r(a) dt Jo 

where T(.) is the Gamma function. From the definition, it becomes apparent that subdiffu- 
sion corresponds to a slowly decaying memory integral in the dynamical equation for p(x, t). 
The operator g^r, < < 2, stands for the Riesz fractional derivative*^ with the Fourier 
transform F{-^pf{x)} = -\k\^f(k). 

The occurrences of the operator oDl~ a and are induced by the heavy-tailed waiting 
times between successive jumps and the heavy-tailed distributions of the jumps, respectively, 



in the underlying continuous-time random walk scheme. The case of a = 1, \x = 2 corre- 
sponds to the standard Fokker-Planck equation. V(x) is an asymmetric periodic potential 



V{x) 



2tt 



sin(27rx) + — sin(47rx) 



(3) 



where Vq and A are the amplitude and the asymmetric parameter of the potential, respec- 
tively. 

Because it is very difficult to solve Eq. (1) analytically and numerically, we used the 
subordinated Langevin method proposed by Magdziarz and coworkers^ to investigate the 
transport. In their method, the solution p(x, t) of Eq. (1) is equal to the probability density 
function of the subordinated process 



Y(t) = X(S t ), 



(4) 



where the parent process X(t) is defined as the solution the stochastic differential equation 



dX(r) 



V'{X{t)) 



(5) 



where L^ir) is the symmetric //-stable Levy motion with the Fourier transform J r {L At (r)} 
e~ T l fc r. Employing the Euler scheme to Eq. (5), one can obtain 



x{ n ) = x{tU, 



X{t ) = 0, 

v 



At + (DAtY^C 



(6) 



(7) 



where % = 1,2, 3... and & are the random variables with standard symmetric //-stable distri- 
bution. The procedure of generating realizations & is the following^ 1 ^ 



sin(//V) 

(cos vy/f 
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where the random variable V is uniformly distributed on(— 7r/2, 7r/2), W has exponential 
distribution with mean one. 

The inverse-time a-stable subordinator S t , which is assumed to be independent of X(t), 
is defined as 

S t = inf{r : U(r) > t}, (9) 

where U{r) is the strictly increasing a-stable Levy motion with Laplace transform 
C{U(t)} = e~ Tka . 

Using the standard method of summing increments of the process U(t) one can get 

[/(r ) = 0, (10) 

U(t j ) = U(t j _ 1 ) + At 1 /%, (11) 
where j = 1, 2, 3... and Q are the skewed positive a-stable random variables^ 1 ^ 5 -. The method 



to generate the random variables is 







sm(a(V + f)) 



cos{V - a{V + § )) 
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(12) 



W 

where V and W have the same definitions as that in Eq. (8). From the above procedures, 
one can obtain the subordinated process Y(t) and its probability distribution function is 
equal to the solution of Eq. (1). For more detailed information on the algorithm, please see 
the Ref. (22). 

In the classical ratchets, one can use the average velocity and effective diffusion coefficient 
to describe the transport. However, for the noise with distribution of a Levy-stable law, the 
mean of the noise and the second moment may do not exist. As a consequence, the classical 
stochastic theory (average velocity and effective diffusion coefficient), which is based on the 
ordinary central limit theorem, is no longer valid. To overcome this problem, Dybiec and 



coworkers^ recently proposed a different approach to the Levy ratchet problem based on 
the quantile line analysis for < /i < 2. 

Quantile line is a very useful tool for investigation of the overall motion of the probability 
density of finding a particle in the vicinity of K^) 18 ' 23 . A median line for a stochastic process 
Y{t) is a function of qo.^it) given by the relationship Pr(Y(t) < qo.sit)) = 0.5. Therefore, 
one can use the derivative of the median to define the group velocity of the particle packet^, 

V a = (13) 

and this definition is valid even for the case of lacking average velocity. 

III. NUMERICAL RESULTS AND DISCUSSION 

In order to investigate the competition between Levy flights and subdiffusion in a ratchet 
potential, we carried out extensively numerical simulations based on the subordinated 
Langevin method 2 ^. For simplicity we set rj = 1.0 and 1 < fi < 2 throughout the work. 
In our simulations, we have considered more than 10 5 realizations to obtain the accurate 
median. In order to provide the requested accuracy of the system the dynamics time step 
was chosen to be smaller than 10~ 3 . We have checked that these are sufficient for the system 
to obtain consistent results. 

Firstly, we will investigate the diffusive properties of the Brownian particles. Usually, 
the types of the diffusion processes are determined by the spread of the distance traveled by 
a random walker. The diffusion is characterized through the power law form of the mean- 
square displacement (x 2 (t)) oc t s . According to the value of the index 5, one can distinguish 
subdiffusion (0 < 5 < 1), normal diffusion (5 = 1) and super diffusion (5 > 1). Here, 
we use the median of square displacement M(x 2 ), instead of mean-square displacement, to 
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characterize the diffusion. Fig. 1 (a) shows the time dependence of M(x 2 )/t for different 
combinations of fi and a without any external potential. It is found that the linear time 
dependence of the median of square displacement, M(x 2 ) oc t, will occur for the case of — =1, 
which indicates the normal diffusion. However, this is not true, for example /x = 1.8 and 
a = 0.9, the process is still non-Markov and non-Gaussian. This pseudo-normal diffusion is 
due to the competition between Levy flights and subdiffusion. Dybiec and coworkers^ have 
presented discussions in detail on this paradoxical diffusion. Fig. 1 (b) presents the time 
dependence of M(x 2 )/t in the presence of a ratchet potential. Interestingly, the pseudo- 
normal diffusion for fx = 1.8 and a = 0.9 with external potentials is not normal. 

Next, we will study the rectified mechanism of the Levy ratchets. Usually, the ratchet 
mechanism demands three key ingredients^ which are (a) nonlinear periodic potential: it 
is necessary since the system will produce a zero mean output from zero-mean input in a 
linear system; (b) asymmetry of the potential, it can violate the symmetry of the response; 
(c)fluctuating: Levy flights can break thermodynamical equilibrium. In Fig. 2 (a), we 
studied the time dependence of the median for different values of the asymmetry parameter 
A at fi = 1.5 and a = 1.0. The median is positive for A > 0, zero at A = 0, and negative 
for A < 0. Therefore, the asymmetry of the potential will determine the direction of the 
transport and no directed transport occurs in a symmetric potential. Now we will give the 
physical interpretation of the directed transport for the case of A = 1. Firstly, the particles 
stay in the minima of the potential awaiting large noise pulse to be catapulted out. The 
particles will be thrown out to the left and the right with the equal probabilities. In this 
case, the distance from minima to maxima is shorter from the right side than that from the 
left side. Consequently, most of the particles are thrown out from the right side, resulting 
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in positive transport. This gives rise to the overall preferred motion to the right. 

Figure 2 (b) gives the time dependence of the median for different combinations of fi and 
a. We find that Levy flights are necessary to obtain the directed transport. For Gaussian 
case (/i = 2.0), directed transport disappears. This is due to the nonthermal character of 
the Levy flights that can break thermodynamical equilibrium. From Fig. 2(a) and (b) we 
can see that the asymmetry of the potential and the non-equilibrium character of the Levy 
flights are the two necessary conditions for directed transport. The direction of the transport 
is determined by the direction of the steeper slope of of the potential and the Levy flights 
can break thermaldynamical equilibrium. These two key ingredients can realize the ratchet 
effects. 

Figure 3 illustrates the dependence of the group velocity V g on the subdiffusive index 
a for different values of the Levy index /i. One can see that group velocity V g increases 
monotonically with the subdiffusive index a. For small values of a, the waiting time between 
successive jumps is very long and it is not easy for particles to pass across the barrier. Thus, 
most particles will stay in their original minima of the potential and the group velocity 
becomes very small. Especially, we also find that there exists a threshold value of a below 
which no directed transport can be obtained. The subdiffusion dominates the transport for 
small values of a (a < 0.7), while the effects of the Levy flights become preponderant for 
large values of a. 

Figure 4 shows the dependence of the group velocity V g on the Levy index ji for different 
values of a. When n — > 2.0, the system is under thermodynamical equilibrium and no 
directed transport appears. For small values of /x, Levy flights are longer and the outliers in 
the Levy noise are larger. In this case, the effects of the asymmetry of the potential become 
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very small, resulting in small group velocity. Therefore, there exists a optimal value of /x at 
which the group velocity takes its maxima. This can also be confirmed by Fig. 3. For very 
small values of a, for example a = 0.5, the group velocity is zero for all values of /x and the 
transport is absolutely dominated by subdiffusion. 

The group velocity V g as a function of noise intensity D is shown in Fig. 5 for different 
combinations of /x and a. The curve is observed to be bell shaped which shows the feature of 
resonance. When D — > 0, the particles cannot pass across the barrier and there is no directed 
current. When D — > oo so that the noise is very large, the effect of the potential disappears 
and the group velocity tends to zero, also. There is an optimal value of D at which the group 
velocity is maximal. There are two intersections (_D c i and D c2 ) between the line of /x — 1.9 
and a = 1.0 and the line of /x = 1.5 and a = 0.9. For simplicity, we define V g (fi,a) as the 
group velocity for different values of /x and a. When D < D cl , V^(1.9, 1.0) > 1^(1.5,0.9), 
Levy flights dominates the transport. When D cl < D < D c2 , 1^(1.9,1.0) < 1^(1.5,0.9), 
the transport is governed by subdiffusion. For the case of D > D c2 , all particles can 
easily pass across the barrier and Levy flights will mainly contribute to the transport and 
1^(1.9,1.0) > 1^(1.5,0.9) 

In Fig. 6, we plot the dependence of the group velocity V g on the amplitude Vq of the 
potential for different combinations of /i and a. When V — > 0, the effects of the potential 
disappear and the group velocity tends to zero. When l/ — > oo, the particles cannot pass 
across the barrier and the group velocity goes to zero, also. Thus, the curve shows a peak. 
Remarkably, there is an intersection between the line of /x = 1.5 and a = 0.9 and the line of 
/x = 1.9 and a = 1.0. This is due to the competition between Levy flights and subdiffusion. 
When V < V c , Levy flights are predominant and 1^(1.5,0.9) > V g (1.9, 1.0). When V > V c , 
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subdiffusion dominates the transport and V g (1.5, 0.9) < V g (1.9, 1.0). In this case, the height 
of the barrier is very high and few particles driven by Levy flights can pass across the barrier, 
the effects of the Levy flights will disappear and subdiffusion will play a major role. 

IV. CONCLUDING REMARKS 

In this paper, we have investigated the directed transport of the Brownian particles in 
a ratchet potential driven by Levy flights coexisting with subdiffusion. We used recently 
developed framework of Monte Carlo simulation 28 which is equal to the solution of the 
fractional Fokker-Planck equation. The group velocity proposed by Dybiec and coworkers 18 
is used to measure the transport. It is found that the group velocity increases monotonically 
with the subdiffusive index, while the group velocity as a function of the Levy index is 
nonmonotonic. The former is caused by the increase of the waiting time between successive 
jumps and the latter is owing to the interplay between Levy flights and the height of the 
barriers. There is a threshold value of a below which the transport is absolutely dominated 
by subdiffusion and the directed transport disappears. The dependences of the group velocity 
on the noise intensity and the amplitude of the potential are also investigated. There is an 
optimal value of the noise intensity (the amplitude of the potential) at which the group 
velocity is maximal. The competition between Levy fights and subdiffusion in the ratchet 
potential is observed on the level of the group velocity as well as the median of square 
displacement. The nonthermal character of the Levy flights and the asymmetry of the 
potential are the necessary conditions for directed transport when the system is in the 
absence of any external driving forces. Because of this competition, we also found the 
pseudo-normal diffusion reported by Dybiec and coworkers^, in which time dependence of 
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the median of square displacement is linear, M(x 2 ) oc t, while the process is still non-Markov 
and non-Gaussian. 

Anomalous transport is becoming widely recognized in a variety of the fields. Beyond 
its intrinsic theoretical interest, the results we have presented may have wide applications 
in some complex systems, such as diffusive transport in plasmas, particles separation with 
non-Gaussian diffusion, and ratchet transport in biology systems that are intrinsically out 
of equilibrium. 
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V. CAPTION LIST 

Fig. I. Time dependence of M(x 2 )/t for different combinations of /x and a: (a)without 
external potential at D = 0.4, solid lines present t * scaling; (b)with external potential 
at D = 0.4, V Q = 5.0, and A = 1.0. 

Fig. 2. Time dependence of the median: (a) for different values of the asym- 
metry parameter A at D = 0.4, V = 5.0, /i = 1.5, and a = 1.0, the inset shows the 
potential profile; (b)for different combinations of /i and a at D = 0.4, = 5.0, and A = 1.0. 

Fig. 3. Group velocity V g versus subdiffusive index a for different values of /i at D — 0.4, 
V = 5.0, and A = 1.0. 

Fig. 4. Group velocity V g versus Levy index /i for different values of a at D = 0.4, 
V = 5.0, and A = 1.0. 

Fig. 5. Group velocity V g as a function of noise intensity D for different combinations of 
ji and a at V = 5.0 and A = 1.0. 

Fig. 6. Group velocity V g versus the amplitude V of the potential for different 
combinations of /i and a at D = 0.4 and A = 1.0. 
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